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Abstract 


A  theory  of  local  and  global  mechanisms  of  heat  generation  and  distribution  in 
carbon-fiber-based  composites  subjected  to  an  alternating  magnetic  field  has  been  proposed.  A 
model  that  predicts  the  strength  and  distribution  of  thermal  generation  through  the  thickness  of 
carbon-fiber-based  laminated  composites  has  been  developed.  Earlier  work  has  established  the 
distribution  of  point  voltages  in  the  plane  of  the  laminate  that  exist  in  the  form  of  potential 
differences  between  fibers  in  adjacent  plies  in  a  cross-ply  or  angle-ply  laminate  system.  In  this 
work,  a  capacitive-layer  microstructure  that  models  the  actual  fiber-reinforced-polymer 
microstmcture  from  a  square-packing  assumption  to  a  series  of  conductive  parallel  plates  is 
formulated.  An  effective  parameter  of  heating,  gamma,  that  establishes  the  distribution  of 
heating  through  the  thickness  is  defined.  Extreme  gradients  in  this  thermal  source  can  exist  with 
peaks  occurring  at  the  interfaces  of  ply-ply  orientation  changes.  An  optimization  study 
establishes  the  effect  of  various  microstmctural  and  macrostractural  parameters  on  the  heating 
parameter,  gamma.  Several  parametric  studies  are  performed  on  the  computer  algorithm,  which 
calculates  gamma  to  further  analyze  these  effects. 
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1.  Introduction 


Induction  heating  involves  the  application  of  an  alternating  magnetic  field  that  subjects  an  object 
to  a  distribution  of  magnetic  flux.  Experimental  observations  have  shown  that  unidirectional  carbon- 
fiber  laminates  subject  to  a  uniform  field  do  not  heat.  In  contrast,  cross-ply  laminates  exhibit 
significant  heating.  Two  mechanisms  of  heating  have  been  identified  that  result  from  an  alternating 
rotational  emf  induced  by  the  magnetic  field:  (1)  joule  losses  in  the  carbon  fibers  induced  by  the 
eddy  currents  and  (2)  dielectric  losses  in  the  polymer.  For  joule  losses  to  dominate  heating,  carbon 
fibers  must  be  in  good  electrical  contact  with  each  other  in  adjacent  layers  of  the  cross-ply  layer. 
Experiments  have  been  conducted  on  cross-ply  laminates  that  incorporate  additional  polymer 
interlayers.  In  this  case,  significant  heating  occurs  for  microstractures  that  eliminate  fiber-fiber 
contact.  This  observation  can  be  explained  by  proposing  that  the  dominate  mechanisms  for  localized 
heating  result  from  dielectric  losses  in  the  matrix  region  separating  fiber-fiber  junctions,  as 
illustrated  in  Figure  1  [1].  In  this  situation,  a  magnetic  field  B  induces  electric  fields  within  the 
carbon  fibers,  creating  a  potential  difference  V  across  the  polymeric  region  of  thickness  h.  A  “planar 
grid”  model  was  developed  to  explore  the  consequence  of  local  heating  resulting  from  dielectric 
losses  in  the  matrix  [2].  hi  this  model,  two  crossing  layers  of  filaments  separated  by  a  polymer  layer 
of  thickness  h  were  treated  as  a  planar  grid  of  capacitors.  The  transfer  of  energy  from  an  applied 
magnetic  field  generated  by  a  Helmholz  coil  was  related  to  the  resulting  distribution  of  potential 
differences  in  the  polymer  plane  separating  the  two  layers.  A  nondimensional  effective  parameter  of 
heating  was  introduced  as  the  voltage  per  unit  magnetic  flux  per  unit  frequency  of  the  applied 
magnetic  field. 

Qualitative  temperature  distributions  over  the  plane  were  observed  by  placing  liquid  crystals  on 
the  surface  of  a  thin-section  cross-plied  carbon/polyetheretherketone  laminate  [2].  These 
observations  showed  local  heating  at  the  fiber  crossover  points,  as  predicted  by  the  planar-grid 
capacitance  model.  Furthermore,  the  temperature  distribution  induced  by  a  circular  coil  was  in 
qualitative  agreement  with  the  predictions.  Laminates  were  fabricated  with  various  thicknesses  of 
polymer  layers  between  the  laminae.  The  local  heating  at  the  fiber-fiber  junctions  persisted. 
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Figure  1.  Induced  Current  Due  to  a  Transverse  (Normal  to  the  Plane  of  Fibers)  Magnetic 
Field  B  in  a  [0/90]s  Cross-Ply.  A  Matrix-Rich  Region  of  Thickness  h  Exists  Between 
the  Orthogonal  Plies  in  Such  a  Laminate.  A  Potential  Difference  V  Exists  Between 
the  Intersecting  Fibers. 


even  though  no  fiber-fiber  contacts  were  possible;  however,  as  the  thickness  of  the  polymer  layer 
was  increased,  the  intensity  of  heating  was  diminished. 

In  this  work,  the  previous  treatment  is  extended  to  account  for  multilayer  interactions  that  occur 
through  the  thickness  in  a  carbon-based  composite.  Previous  work  modeled  the  interaction  between 
single  transverse  fibers  using  an  electric  field  concept  [  1  ]  and  the  planar  distribution  of  current  based 
on  an  equivalent-resistance  assumption  [2]  resulting  in  a  qualitative  view  of  thermal  generation.  In 
the  present  treatment,  it  is  shown  that  the  multilayer  systems  can  give  rise  to  complicated  potential 
differences  between  the  fibers  through  the  thickness  of  the  laminate.  First,  a  capacitive  layer  analogy 
is  developed  and  simplifications  are  introduced  to  provide  estimates  of  the  effective  fields  within  the 
polymeric  region  that  can  be  consolidated  into  an  effective  thickness  parameter  for  the  polymeric 
material.  Through  this  approach,  the  effects  of  a  complex  distribution  of  potential  differences  can  be 
captured  in  the  form  of  a  nondimensional  heating  parameter.  Parametric  studies  based  on  this  model 
identify  the  role  of  various  microstructural  features  in  the  induction  welding  of  carbon-fiber-based 
thermoplastic  composites.  Experimental  evidence  is  presented  to  support  these  predictions. 
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2.  Formulation 


Figure  2  illustrates  an  idealized  system  of  two  orthogonal  plies,  each  containing  three  fiber 
“layers,”  separated  by  a  distance  ho  =  2hi,  as  shown.  In  accordance  with  the  planar  grid  model  [2], 
each  fiber  layer  pair  between  the  adjacent  plies  of  the  2-ply  system  is  initially  assumed  to  have  the 
same  potential  difference.  For  example,  the  potential  difference  between  fiber  layer  no.  3  in  the  0® 
ply  and  fiber  layer  no.  3  in  the  90°  ply  is  initially  assumed  to  be  equivalent  to  the  potential  difference 
between  fiber  layer  no.  3  in  the  0°  ply  and  fiber  layer  no.  2  in  the  90°  ply.  Each  of  the  n^  =  3^  =  9 
fiber  layer  interactions  is  initially  assumed  to  have  the  same  potential  difference  as  calculated  by  the 
planar-grid  model.  This  implies  that  each  fiber-fiber  interaction  creates  a  complicated  electric  field 
between  two  nonparallel  conducting  cylinders.  As  an  example,  for  a  20-cm-square  (~8  in.)  2-ply 
cross-ply  laminated  plate,  there  are  approximately  2.5e  +11  interactions  to  consider.  Consequently, 
a  simplified  model  is  needed  to  estimate  the  electric  field  in  the  polymeric  regions  and,  in  doing  so, 
define  an  effective  parameter  for  the  thickness  of  polymeric  material,  h,  through  which  the  field  acts. 
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Figure  2.  Simplified  Representation  of  the  Cross  Section  of  a  [0/90]  Cross-Ply  Laminate. 

Electric  field  lines  between  fibers  act  almost  entirely  in  the  through-thickness  sense,  interacting 
with  nearly  all  of  the  polymeric  volume  in  the  region  between  the  two  fibers.  Under  these 
conditions,  it  is  plausible  to  represent  rows  of  fibers  as  contiguous  conductive  material. 
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Furthermore,  the  fibers  can  be  taken  to  be  perfect  conductors  in  comparison  to  the  insulating 
polymer.  With  the  assumption  that  the  fiber  volume  fraction  and  laminate  thickness  are  continuous 
in  an  elemental  planar  region,  the  electromagnetic  interaction  between  two  arbitrary  fibrous  layers 
can  be  treated  as  conductive  plates  separated  by  some  particular  distance,  h. 

Initially,  consider  just  two  such  interacting  conductive  plates  making  up  a  “capacitive  layer”  pair 
in  the  x-y  plane.  A  potential  difference  exists  between  the  plates  that  is  a  function  of  the  planar 
coordinates,  x  and  y: 


V(x,y)  =  -jE(x,y)-d^  =  E(x,y)h,  (1) 

where  E(x,y)  is  the  electric  field  acting  through  the  distance,  h,  separating  the  two  plates.  Therefore, 
a  charge  distribution  exists: 


q(x,y)  =  e„KjE(x,y)-dS,  (2) 

where  q(x,y)  is  the  coulombic  charge,  Eo  is  the  permittivity  of  vacuum,  k  is  the  relative  dielectric 
constant  of  the  material,  and  S  is  the  Gaussian  surface  vector.  If  some  elemental  surface  area  AA  is 
considered, 


q(x,y)A  =eoK<l>E  =eoKE(x,y)AA 


(3) 


or 


q(x«y)AA  = 


£„KAAV(x,y) 

h 


(4) 


If  AA  is  taken  large  enough  such  that  its  linear  dimensions  are  much  greater  than  h,  the  capacitance 
at  the  point  (x,y)  can  be  obtained: 
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(5) 


r/'v  v'k  _  q(x,y)  _  eoKE(x,y)AA  _  e^kAA 
^  V(x,y)  E(x,y)h  h 

The  heating  at  point  (x,y)  can  now  be  written  as 

W(x,  y)^^  =  CO  tan  6  C(x,  y)^  V(x,  y)^ ,  (6) 


where  (o  is  the  angular  frequency  of  the  electric  field  and  tan5  is  the  imaginary  part  of  the  complex 
dielectric  constant  for  the  polymer.  Incorporating  equation  (5)  for  the  capacitance  yields 


W(x,y)^  =  (0tan9-^^^^  V(x,y)^ 
n 


pV(x,y)^ 

h 


(7) 


Note  that  the  quantity  p/h  is  equivalent  in  units  to  1/R,  where  R  is  a  resistance  indicating  that 
dielectric  heating  follows  a  joule-type  loss  relationship.  In  the  planar  grid  model  of  Fink, 
McCullough,  and  Gillespie  [2],  fairly  large  (2.5-cm  square)  element  sizes  (AA)  were  determined  to  be 
satisfactory  representatives  of  the  potential  distribution  for  centered-coil  problems.  It  was  shown  that 
the  gradient  in  V(x,y)  in  the  plane  is  much  steeper  for  comer-  or  edge-placed  coil  arrangements  due 
to  the  nature  of  eddy  current  formation  in  cross-ply  laminates.  Comer-  or  edge-placed  coil 
arrangements  require  a  finer  mesh  for  accurate  in-plane  voltage  analyses. 

As  noted  previously,  a  through-thickness  model  must  account  for  many  paraUel-plate  capacitive 
interactions.  Figure  3(a)  illustrates  three  planar  fiber  layers  above  the  0790°  interface  and  three 
below.  The  matrix  interface  region  between  the  0  and  90  is  assigned  a  thickness  of  ho;  the  thickness 
of  the  matrix  layer  within  both  the  0  and  90  regions  is  assigned  the  thickness  d*.  The  quantity  d*  is 
dependent  upon  the  volume  fraction  of  fibers  in  a  ply  and  the  packing  geometry.  For  this  example, 
nine  pair-wise  fiber  layer  interaction  combinations  are  possible.  Each  isolated  pair  initially  carries  the 
same  potential  difference,  V(x,y);  however,  the  separation  distance,  h,  varies  with  each.  For  each 
isolated  pair,  the  equation  for  heating,  equation  (6),  holds.  However,  as  the  isolated  pairs  are  brought 
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Figure  3.  Series  of  Capacitive  Layers  Representing  (a)  a  [0/90]  Cross-Ply  Laminate  in  Which 
Each  Ply  Consists  of  Only  Three  Fibers  Through  the  Thickness  on  Average  and  (b) 
a  Parallel  Plate  Capacitor  With  a  Conductor  Replacing  Some  of  the  Dielectric 
Material. 


together  to  form  the  overall  layered  structure,  the  capacitance  and  voltage  effectively  change  due  to 
the  presence  of,  N,  intervening  conducting  fiber  layers,  each  of  thickness  df*,  between  the  fiber 
layer  pair  under  consideration. 

The  introduction  of  a  conductive  slab  into  a  capacitor  effectively  decreases  the  distance  between 
the  charged  plates.  Figure  3(b)  illustrates  the  situation  in  which  a  conductor  is  introduced  within  a 
capacitor  containing  a  dielectric.  A  is  the  area  of  the  charged  plates  in  the  capacitor,  and  q  is  the 
charge  on  each  plate.  Before  the  conductive  slab  is  introduced. 
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After  the  slab  is  introduced, 


C 

new 


EqAk. 


'dd 


r.  bnV 


1— 


(10) 


and 


(11) 


where  bo  is  defined  in  Figure  3(b)  as  the  thickness  of  a  conductive  slab  that  replaces  dielectric  material 
in  a  parallel-plate  capacitor.  The  capacitance  is  increased  by  the  same  amount  that  the  potential 
difference  is  decreased.  However,  since  heating  is  proportional  to  CV^,  the  heating  rate,  from 
equation  (6),  is  decreased  by  a  factor  of  (1  -  bo/h).  This  decrease  in  heating  can  be  accounted  for  by 
an  effective  (1/h)  value  for  each  possible  fiber-fiber  interaction  through  the  thickness.  In  the 
following  discussion,  the  quantity  bo  of  Figure  3(b)  is  equivalent  to  the  summation  of  effective 
thicknesses  of  fiber  layers  (conductive  slabs)  separating  the  two  fiber  layers  under  consideration. 


The  effect  of  several  interactions  can  be  illustrated  for  the  [0/90]  laminate  of  Figure  3(a).  The 
plies  of  this  laminate  are  assumed  to  have  only  three  fiber  layers  for  this  discussion.  The  “upper”  fiber 
layers  are  labeled  Ui,  U2,  and  ua  and  the  “lower”  fiber  layers  are  labeled  h,  I2,  and  I3.  The  thickness  of 


the  region  between  the  0°  and  90°  plies  is  labeled  ho.  The  polymeric  regions  are  labeled  j  =  1, 2,  or  3. 
Table  1  summarizes  the  nine  possible  fiber-fiber  interactions  and  the  values  of  1/h  for  each  polymeric 
region  between  fiber  planes.  The  quantity  d*  is  the  capacitive  layer  model’s  value  for  the  distance 
between  conductive  plates  (i.e.,  the  thickness  of  the  polymer  regions). 


Table  1.  Individual  Inverse  Fiber-Fiber  Interaction  Distances  for  an  Ideal  Cross-Ply 
Laminate  in  Which  Each  Ply  Consists  of  Three  Fiber  Layers  Through  tiie  Thickness 


(1.1) 

(1.2) 

(1.3) 

(2,1) 

(2,2) 

(2,3) 

(3.1) 

(3,2) 

(3.3) 

Total  Per  Polymeric 
Region 

U2  to  U3 

j  =  3 

d*/B 

d*/C 

d*/D 

2  2+a  H  * 

11/ .  .^2 

‘=2i=«(ho-i-id  ) 

Ui  to  U2 

j  =  2 

d*/A 

d*/B 

d*/C 

d*/B 

d*/C 

d*/D 

Ui  to  li 

j=l 

1/ho 

ho/A 

ho/B 

ho/A 

hc/B 

ho/C 

hc/B 

ho/C 

ho/D 

li  to  I2 
j  =  2 

d*/A 

d*/B 

d*/B 

d*/C 

d*/C 

d*/D 

I2  to  I3 

j  =  3 

d*/B 

d*/C 

d*/D 

2  2+a  (J* 

-=2i-(ho+id  ) 

2  2+a  1 

Total:  a5S(ho+id*) 

Notes;  A=(ho  +  d*)^ 

B  =  (ho  +  2d*)^ 

C  =  (ho  +  3d*)^aDd 
D=(ho  +  4d*)l 


The  distance  between  ui  and  li  is  ho  and,  since  there  are  no  other  fibrous  conductive  layers 
between  them,  1/h  =  1/ho.  For  the  interaction  between  Ui  and  I2,  labeled  as  ( 1 ,2)  in  the  table,  there  are 
two  polymeric  regions  to  consider  but  only  one  fiber-fiber  interaction.  The  total  distance  between  Ui 
and  I2  is  ho -f  d*  so  the  total  1/h  for  this  interaction  is  l/(ho  -1-  d*).  However,  the  heat  generation  is 
divided  between  two  polymeric  regions;  one  region  has  thickness  ho  and  the  other  has  thickness  d*. 
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Therefore,  the  fraction  of  heating  in  each  is  represented  as  ho/(ho  +  d*)  and  d*/(ho + d*),  respectively. 
The  total  heating  in  each  region,  then,  is  given  by 


and 


^_ho _ i_ 

ho+d*ho+d* 


d*  1  _  d* 

^h>j=2  ho+d  ho+d  (ho+d*y 


(12) 


(13) 


for  the  j  =  1  and  j  =  2  regions,  respectively,  of  the  (1,2)  interaction.  The  remainder  of  Table  1  is 
completed  accordingly.  The  column  labeled  ‘Total  Per  Polymeric  Region”  gives  the  total  effective 
inverse  distance  for  each  polymer  region  taking  into  account  all  fiber-fiber  interactions.  The 
summation  of  these  terms  gives  the  total  effective  inverse  h  for  the  ply-ply  interaction.  From  here  on, 
this  summation  is  termed  the  “gamma  parameter”  and  is  given  the  symbol  yj  for  the  j'*  polymeric 
region  through  the  thickness. 

This  example  gives  the  results  for  a  laminate  with  only  three  fiber  layers  above  and  below  the  ply- 
ply  interface.  The  general  case  for  a  two-ply  laminate  consisting  of  m  fiber  layers  in  the  “top”  ply  and 
n  fiber  layers  in  the  “bottom”  ply  is  given  in  equations  (14)  and  (15)  for  the  effective  inverse  h 
parameters  yj  and  y,otoi,  respectively. 


fi  Y'" 

h  >j=i 


m-ln+a-l 

=  11 


(ho+id*)' 


-,  and 


(14) 


^  /j{n') 


1  A  m-1  n+a-l  d 

V  =11  7 - 

n  jj(m)  »=H  i=d  (ho  +  id  j 


(14a) 


where  j  =  1, ...,  m; 
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(14b) 


fi 

^  ^  Jj(n) 


a  ^ 

h* 


n-1  m+a-1  q 

~  S  S  ~/  ,  V2  ’ 

;j(n)  ‘‘=H  (ho+id  / 


where  j  =  1, ...,  n;  and 


eff 


y,  ^  /total 


a  ^ 

h' 


/total 


m-lii+a-1  1 

=  II  TT — 

a=0  i=:a  ^ho+ld  ) 


(15) 


Equations  (14a)  and  (14b)  are  identical  if  m  =  n. 

Quantities  of  special  interest  are  the  thermal  energy  (heating)  at  the  interface  of  two  adjacent 
plies  represented  by  yj ,  the  total  heating  in  the  two-ply  laminate  represented  by  7totai>  and  the  ratio  of 
these  values,  which  represents  the  fraction  of  total  heating  that  occurs  in  the  interface  region.  To 
qualitatively  evaluate  the  role  of  these  parameters,  it  is  useful  to  replace  equations  (14)  and  (15)  by 
integral  forms;  viz., 


J  (m-l)(n+a-l) 

U  a=0  i=a  + 1  j 


dida. 


where  b  =  ho/d*.  Equation  (16)  can  be  approximated  as 


(b -I- m  -  iXb  H- n  - 1) 
b(b  -h  n  -I-  m  -  2) 


For  the  case  where  m  =  n,  this  can  be  further  reduced  to 


7i  = 


_1_,  1  (b  +  n-lX 
d*  b(b-h2n-2) 


(16) 


(17) 


(18) 


Likewise,  for  the  total  effective  inverse  h. 
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(19) 


7,otai  =  ^[(rlnr  -  r)-  (sins  -  s)-  (tint  - 1)+  (blnb  -  b)], 
d 

where 

r  =  b  +  n  +  m-2, 
s  =  b  +  m  -  1,  and 
t  =  b  +  n  -  1. 


For  the  case  where  m  =  n, 

Ttotai  =:^(gl^ig-2hlnh  +  blnb),  (20) 

d 

g  =  b  +  2n  -  2,  and 
h  =  b  +  n  -  1. 

The  optimization  of  overall  energy  production  is  obtained  by  the  maximization  of  equations  (17)  and 
(19).  These  conditions  yield  the  requirement  that  m  =  n  be  maximized  and  ho  be  minimized. 


If  equations  (18)  and  (20)  are  divided  by  the  volumes  over  which  they  act  and  the  ratio  of  the  two 
taken,  the  ratio  of  heating  is  obtained  at  the  interface  to  heating  in  the  laminate  (equation  [21]).  In 
welding  operations,  it  may  be  desirable  to  preferentially  (or  nonpreferentially)  heat  the  interface  so 
this  ratio  should  be  as  large  (or  small)  as  possible. 


7i 

7  total 


b 


2Z-- 


7i 


b  +  2n 


(21) 


11 


where 


Z  =  — ^ —  (n  -  l)ln 
b  +  2n^  ^ 


b  +  2n-2^ 

b+n-1  J 


Maximization  requires  that  ply  thicknesses,  m  =  n,  be  maximized  and  that  the  interface  thickness,  ho, 
be  minimized  as  before.  Therefore,  increasing  heating  by  either  increasing  ply  thickness  or 
decreasing  ho  also  increases  the  preferential  nature  of  heating  at  the  interface. 

As  implied  earlier,  the  quantity  d*  is  related  to  the  volume  fraction  fiber  (Xf)  and  the 
fiber-packing  geometry.  Several  ideal  packing  geometries  exist;  the  most  common  are  rectangular, 
square,  and  hexagonal.  A  square-packing  assumption  was  used  for  this  study  for  simplicity  and 
adaptability  to  the  capacitive-layer  model.  As  a  caveat  for  this  assumption,  rectangular  and 
hexagonal  packing  geometries  were  analyzed  separately  and,  with  deference  to  the  more 
straightforward  square-packing  geometry,  were  found  to  have  no  significant  influence  on  the 
through-thickness  potential  difference  results  when  realistic  thicknesses  of  composite  laminates  were 
considered.  Figure  4(a)  shows  the  parameters  involved  and  Figure  4(b)  shows  the  modeled  system 
with  its  corresponding  parameters.  The  quantity  a  is  the  dimension  of  the  unit  cell  and  is  equivalent 
to  the  summation  of  the  fiber  diameter  df  and  the  distance  between  fibers  d.  The  quantities  df*  and 
d*  are  the  fiber  diameter  and  fiber  separation  dimensions  corresponding  to  the  model’s  unit  cell.  The 
unit  cell  dimension  a  is  taken  to  be  constant,  establishing  the  condition  that 

d+df=d*+d;  (22) 

and  the  fiber  volume  fraction  remains  constant,  requiring  that, 


d; 

d*  H-d;  ■ 


(23) 
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Figure  4.  Unit  Cells  for  (a)  Square-Packing  Assumption  and  (b)  Capacitive-Plate  Assumption. 


Equation  (23)  can  be  used  to  specify  d*  in  terms  of  the  fiber  volume  fraction  Xf.  Combining 
Equations  (22)  and  (23)  gives  d*  in  terms  of  the  microstructural  parameters  d  and  df,  . 

d*  =  (l-XfXd-i-df).  (24) 


For  the  square-packing  geometry. 


_ 

”4(d  +  d,y  ■ 


(25) 


Combining  these  relationships  yields  the  effective  interlayer  distance  d*  in  terms  of  the 
microstructural  parameters: 


2^ 


(26) 


The  question  arises  as  to  how  many  fiber  layers  are  necessary  for  the  square-packed/capacitive-layer 
model  to  be  appropriate.  The  condition  of  uniform  fiber  volume  fraction  (equation  [23])  is  not 
necessarily  met  when  the  number  of  fiber  layers,  n,  is  low.  Equation  (27)  gives  the  ratio  of  effective 
polymer  thickness  to  total  thickness  for  the  square-packed  model: 


leff  _  (n  ~ 
t  (n-l)d’-hnd;  ■ 


(27) 


Figure  5  shows  the  plot  of  equation  (27)  (with  Xf  =  0.60)  vs.  n,  indicating  a  convergence  to  the 
expected  value  of  0.40.  APC-2  has  a  fiber  volume  fraction  of  approximately  0.60  and  each  prepreg 
ply  has  between  16  and  20  fiber  layers  upon  consolidation.  The  results  of  Figure  5  indicate  that  it  is 
plausible  to  use  the  capacitive-layer  model  for  laminates  of  nominally  one  ply  and  greater. 


Number  of  fiber  layers,  n 

Figure  5.  Plot  of  Equation  (27)  Showing  the  Convergence  of  the  Ratio  of  Effective  to  Actual 
Thickness  in  the  Square-Packed  Capacitive-Layer  Model  for  Fiber-Fiber 
Interactions  Through  the  Thickness  of  a  Laminate  Consisting  of  n  Fiber  Layers. 
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3.  Model  Predictions 


Figure  6  shows  an  outline  of  a  computer  algorithm  used  to  determine  the  effective  heating 
parameter  7.  Blocks  1  and  2  are  based  on  the  relationships  previously  discussed.  Block  3  is  the 
summation  of  the  ply-ply  interactions.  Most  case  studies  were  performed  using  laminates  with  only 
one  change  in  ply  orientation  through  the  thickness.  The  program  calculates  yjd*,  where  j  is  defined 
in  Figure  3.  The  quantity  b  =  ho/d*  is  used  so  that  all  calculated  quantities  are  in  terms  of  the 
model’s  capacitive-layer  separation  distance  d*,  and  the  result  is  nondimensional. 


Laminate  Characteristics: 

Material  Characteristics:  •  #  m  layers  above  interface 

•  Fiber  Volume  Fraction  •  #  n  layers  below  interface 

•  Fiber  Diameter  •  interply  matrix  <rich  region 

I  of  thickness,  h^ 


Figure  6.  Outline  of  Capacitive-Layer  Computer  Algorithm. 


3.1  Effect  of  Interply  Thickness.  Three  relevant  cases  for  polymer  thickness  of  ho  =  125  pm, 
25  pm,  and  3.2  pm  are  presented  as  examples  of  model  prediction.  Table  2  provides  the  program 
input  parameters  for  the  three  cases.  The  calculations  described  to  this  point  do  not  take  into  account 
the  distribution  of  heating  per  unit  thickness  of  the  regions  j.  Since  region  j  =  1  is  of  thickness  ho  and 
regions  j  =  2, . . .,  n  or  m  are  of  thickness  d*,  the  normalized  yd*  values  are  used.  The  ratio  of 
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interfacial  heating  to  total  heating  in  the  ply-ply  interface  region  is  given  in  Table  3  for  each  case.  As 
the  interply  matrix-rich  region  decreases,  ytotai  increases,  as  expected,  and  ji  decreases.  Heating, 
however,  is  a  function  of  the  volume  of  material  to  which  energy  is  being  supplied;  so,  as  the  third 
and  fourth  columns  of  Table  3  indicate,  the  normalized  output  representing  heating  increases  for  the 
entire  laminate  (at  even  faster  rates)  with  decreasing  interply  thickness,  ho-  The  fifth  column 
indicates  the  ratio  of  interface  to  total  laminate  heating,  as  defined  in  equation  (21).  The  high  ratio 
drops  off  rapidly  with  increasing  interply  thickness.  Consequently,  maximum  interfacial  heating  is 
expected  to  occur  when  the  thickness  of  the  laminate  is  maximized  and  the  thickness  of  the 
interfacial  polymer  is  minimized. 


Table  2.  Input  Parameters  for  Case  Studies  of  Fiber-Fiber  Submodel.  Xf  =  0.60,  df  =  7  pm, 
20  Fiber  Layers  Per  Ply,  and  d*  =  3.2  pm. 


Case 

0®  Ply  Layers,  m 

90°  Ply  Layers,  n 

Interply  Thickness,  ho 

b  =  ho/d* 

1 

60 

60 

40 

2 

60 

60 

25  pm  (1  mil) 

8 

3 

60 

60 

3.2  pm  (1/8  mil) 

1 

Table  3.  Comparison  of  ytotai  and  yi  for  the  Example  Cases  in  the  Fiber-Fiber  Submodel. 

Normalization  Occurs  Against  the  Ratio  of  Interface  Thickness  to  Effective  Fiber 
Layer  Spacing  From  Table  2. 


Case 

b  =  ho/d* 

(1) 

(2) 

(3) 

(4) 

(5) 

Yid* 

Yiotald* 

Normalized 

7id* 

Normalized 

Ytotald* 

Yl/Ytotal 

Ratio  of  Col. 

(3)  to  (4) 

1 

40 

18.31 

28.66 

0.46 

0.18 

2.55 

2 

8 

12.95 

65.39 

1.62 

0.57 

3.18 

3 

1 

4.98 

82.68 

4.98 

0.68 

7.32 

A  parametric  study  was  performed  to  ascertain  the  effect  of  varying  the  interply  thickness,  ho, 
while  maintaining  a  constant  ply  thickness  of  a  single  ply  on  either  side  of  the  ply-ply  interface 
([0/90]).  Values  of  df  =  7  pm,  =  0.60,  and  m  =  n  =  20  were  input  into  the  model  with  various 
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values  of  b  =  ho/d*.  The  values  for  fiber  diameter  and  fiber  volume  firaction  are  within  the  range  of 
values  reported  by  the  manufacturer  and  in  the  literature  [3-5]  and  the  number  of  fibers  through  the 
thickness  is  within  the  range  of  values  determined  in  a  micrographical  study  [6].  Graphical  results 
are  presented  in  Figures  7  and  8.  Figure  7  plots  nondimensionalized  “energy”  input  to  the  total 
laminate  and  the  “center”  interply  region  vs.  nondimensionalized  interply  thickness.  The  series 
solution  is  the  exact  solution  of  the  model.  The  integral  solution  is  provided  to  qualify  its  use  in  the 
optimization  study. 

Note  that,  at  low  interply  thickness  (b<~'13)  in  Figure  7,  the  energy  production  in  the  interply 
region  increases  with  increasing  ho  then  decreases  slowly  past  that  point. 

From  equation  (21),  the  ratio  of  total  energy  input  to  the  resin-rich  interply  region  to  total  energy 
input  to  the  entire  laminate  was  examined  and  found  to  increase  as  the  interply  region  increased. 
Figure  7  indicates  this  trend  since  the  total  energy  input  into  the  interply  region  is  significantly  lower 
than  that  of  the  laminate  and  the  ratio  decreases  significantly  with  increasing  interply  thickness.  This 
difference  is  shown  quantitatively  in  Figure  8  as  the  ratio  of  specific  energy  or  heating  in  the  interply 
region  to  that  in  the  laminate  (per  unit  thickness).  Figure  8  illustrates  the  important  relationship  of 
increasing  localization  of  heating  with  decreasing  interply  resin  thickness. 

3.2  Varying  Ply  Thickness.  Increasing  the  thickness  of  the  0°  and  90°  plies  within  the  laminate 
has  the  opposite  overall  effect  of  independently  increasing  the  interply  thickness.  The  results  of  a 
parametric  study  involving  several  [On/OOn]  cases  are  graphically  depicted  in  Figures  9  and  1 0,  where 
the  y-axis  variables  are  equivalent  to  Figures  7  and  8,  respectively,  of  the  last  section. 

A  comparison  of  Figure  10  with  Figure  8  shows  that  it  is  beneficial  to  decrease  the  thickness  of 
the  polymer  interply  region  while  increasing  the  number  of  plies  for  maximized  interfacial  heating. 
A  slight  decrease  in  ho  can  possibly  produce  more  significant  increases  in  heating  than  can  a  larger 
magnitude  increase  in  the  ply  thickness.  Figure  8  suggests  that  heating  can  be  controlled  by  careful 
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Figure  7.  Results  of  Parametric  Study  of  Capacitive-Layer  Model  to  Determine  Effects  of 
Increasing  Interply  Resin  Thickness  on  Dimensionless  Heating  Parameter. 
Comparison  of  Interface  Energy  Input  (y/d*)  to  Total  Energy  Input  (ytotaid*)  for 
Both  the  Exact  Series  Solution  and  the  Approximate  Integral  Solution. 


Figure  8.  Results  of  Parametric  Study  of  Capacitive-Layer  Model  to  Determine  Effects  of 
Increasing  Interface  Resin  Thickness  ho  on  Dimensionless  Heating  Parameter  y. 
Ratio  of  Center  to  Average  Interfiber  Heating,  yi/[(yiotai-yi)/(m+n)],  Plotted  Against 
Increasing  Interface  Thickness.  This  Represents  the  Multiplicity  of  Average 
Interfiber  Heating  That  Is  Consumed  by  the  Interply  Resin.  For  Example,  at  b  =  1, 
Which  Represents  an  Interply  Thickness  Equal  to  d*,  the  Heating  in  the  Interply 
Resin  Rich  Region  Is  Approximately  Six  Times  the  Average  Heating  in  the  Rest  of 
the  Laminate. 
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Half-Thickness  in  Number  of  Plies 

Figure  9.  Comparison  of  Interface  “Enei^  Input”  (yid*)  to  Total  “Energy  Input”  (ytotaid*) 
for  Both  Exact  Series  Solution  and  Approximate  Integral  Solution. 


Half-Thickness  in  Plies 

Figure  10.  Ratio  of  Center  to  Average  Interfiber  “Heating,”  yi/[(ytotai-yi)/(m  +  n)],  Plotted 
Against  Increasing  Ply  Thickness.  This  Represents  the  Multiplicity  of  Average 
Interfiber  Heating  That  Is  Consumed  by  the  Interply  Resin.  For  Example,  at  a 
Laminate  Thickness  of  2  Plies  (Half-Thickness  =  1),  the  Heating  in  the  Interply 
Resin-Rich  Region  Is  Approximately  2.3  Times  the  Average  Heating  in  the  Rest  of 
the  Laminate. 
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control  of  the  thickness  of  the  resin-rich  interply  region — a  relatively  easy  parameter  to  control  in 
processing,  although  with  possible  detrimental  effects  on  mechanical  properties  for  large  h^j. 

4.  Discussion 

While  the  model  presented  in  this  paper  determines  the  through-thickness  distribution  of  the 
effective  parameter  of  heating  y,  it  must  be  combined  with  the  planar  voltage  distribution  model  [2], 
the  local  heat-generation  model  [1],  and  a  finite  element  heat-transfer  model  [6],  which,  together, 
yield  the  thermal  distribution  of  heating  at  the  surface  of  the  laminate.  Extensive  testing  that  strongly 
supports  the  combined  models’  predictions  has  been  performed.  The  details  of  the  experimental 
procedure  used  and  a  comprehensive  verification  of  the  global  model  [6]  have  been  provided 
separately.  A  full  discussion  of  the  testing  procedure  and  the  other  models  is  not  within  the  scope  of 
this  paper;  however,  a  few  results  are  provided  to  initially  verify  the  through-thickness  model  as  a 
component  of  the  global  model,  as  well  as  the  validity  of  the  dielectric  loss  mechanism. 

The  two  key  parameters  discussed  previously  are  the  polymeric  interface  thickness  and  the  actual 
ply  thickness  (or  number  of  plies)  in  a  ply-ply  interaction.  The  predicted  effect  of  varying  the 
interface  thickness  on  the  heating  parameter  y  (Figure  8)  indicated  a  rapid  decrease  in  y  with 
increasing  interface  resin  thickness  ho.  Figure  1 1  compares  predicted  and  experimental  steady-state 
surface  temperatures  at  one  point  for  [0/90]  specimens  manufactured  with  varying  thicknesses  of 
interply  polymer.  The  experimental  error  bars  (heavy  lines)  are  indicative  of  the  ±1  °C  reading  of  the 
digital  thermocouple  and  an  additional  ±5%  error  on  the  accuracy  of  the  placement  of  the 
thermocouple.  The  theoretical  results’  error  bars  (light  lines)  are  the  results  of  the  application  of 
upper  and  lower  bounds  on  all  known  parameters,  including  geometrical  and  material  property 
characteristics.  Similar  results  in  consideration  of  the  ply  thickness  are  shown  in  Figure  12.  While 
these  results  do  not  provide  a  direct  proof  of  the  proposed  model,  they  do  provide  persuasive 
evidence  supporting  the  combined  models’  foundation. 


20 


Equilibrium  Terr^erature  (degrees  Celsius) 


Nominal  Interface  Resin  Thickness  (mil) 

Figure  11.  Comparison  of  Predicted  and  Observed  Surface  Temperatures  Showing  Effect  of 
Increaring  Interface  Thickness  (h©)  on  Equilibrium  Surface  Temperature.  The 
Heavy  Lines  Represent  the  Experimental  Results  While  the  Thin  Lines  Represent 
Predicted  Results  With  Upper  and  Lower  Bounds. 
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Equilibrium  Temperature  (°C) 


Number  of  Plies  on  Each  Side  of  Interface  (M  =  N) 

Figure  12.  Comparison  of  Predicted  and  Observed  Surface  Temperatures  Showing  Effect  of 
Increasing  the  Number  of  Plies  on  Each  Side  of  the  Interface.  The  Heavy  Lines 
Represent  the  Experimental  Results  While  the  Thin  Lines  Represent  Predicted 
Results  With  Upper  and  Lower  Bounds. 
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Other  observations  further  support  the  dominate  role  of  dielectric  losses  rather  than  joule  losses  in 
the  induction  heating  of  carbon-fiber  composites.  As  indicated  previously,  observed  heating  patterns 
show  significant  heating  only  where  conductive  paths  overlap.  Different  polymers  were  used  to 
separate  carbon-fiber  laminae.  Heating  characteristic  changed  with  the  change  m  polymer  interlayer. 
This  supports  the  proposal  that  heating  is  due  to  dielectric  losses  in  the  polymer  since  the  dielectric 
behavior  of  the  polymers  (PEEK  and  PEI)  have  different  dielectric  loss  spectra.  An  alternate 
explanation  for  the  heating  behavior  could  be  attributed  to  dielectric  breakdown  in  the  thin  polymer 
layers  between  fibers.  Such  a  breakdown  would  provide  conducting  paths  to  connect  the  conducting 
carbon  fibers.  However,  cycle  tests  show  no  change  in  heating  rate  between  the  first  cycle  and  later 
cycles.  If  dielectric  breakdown  had  occurred,  the  later  cycles  should  have  exhibited  different  heating 
characteristics.  The  fact  that  rapid  heating  is  attainable  at  relatively  large  interply  thicknesses  further 
discounts  the  role  of  dielectric  breakdown,  as  well  as  charge  transfer  mechanisms. 


5.  Summary 


A  capacitive-layer  model  has  been  developed  to  establish  the  mechanisms  for  fiber-layer 
interaction.  The  model  uses  a  capacitive-layer  parallel-plate  assumption  to  develop  the  equations  that 
describe  the  electromagnetic  interaction  between  fibers  in  adjacent  plies  and  the  resulting  reaction  of 
the  interply  and  interfiber  polymeric  regions  to  the  induced  electric  fields.  This  model  effectively 
incorporates  the  key  microstructural  (fiber  diameter,  fiber  volume  fraction)  and  macrostructural 
(interply  thickness,  laminate  thickness)  characteristics  of  the  laminated  plate.  Each  of  these 
parameters  affects  the  heating  in  a  manner  independently  predicted  by  the  capacitive-layer  model. 
The  model  utilizes  a  nondimensional  parameter,  which  represents  the  thermal  generation  within  a  ply- 
ply  interaction.  This  parameter,  yd*,  was  then  normalized  with  respect  to  volume  for  comparison 
with  other  systems  and  between  polymeric  regions  within  the  laminate.  Note  that  the  model  initially 
assumes  that  each  opposing  fiber-fiber  pair  is  separate  fi-om  its  neighbors  and  carries  an  equivalent 
voltage,  Vxy.  When  the  model  superimposes  the  system  of  fibers  into  its  actual  form,  the  actual 
voltage  distribution  is  proportional  to  the  square  root  of  y  times  the  equivalent  voltage  Vxy. 
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Parametric  studies  identified  the  desired  extremes  of  key  parameters  for  overall  thermal  generation 
and  preferential  (or  nonpreferential)  heating  in  the  ply-ply  interface  region.  These  studies  showed 
that,  in  order  to  maximize  the  heating  at  any  point  in  the  plane  of  the  cross-ply  or  angle-ply  laminate, 
the  fiber  volume  fraction  should  be  maximized,  the  fiber  diameter  should  be  minimized,  the  interply 
resin  thickness  should  be  minimized,  and  the  ply  thickness  above  and  below  the  interface  should  be 
maximized.  This  model  has  been  indirectly  verified  through  a  global  model  analysis  and  the 
execution  of  a  rigorous  experimental  matrix. 
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J  WDLIS 

3815  MEDFORD  ST 

LOS  ANGELES  CA  90063-1900 

1  ADVANCED  GLASS  FIBER  YARNS 

T  COLLINS 

281  SPRING  RUN  LN  STE  A 
DOWNINGTON  PA  19335 

1  COMPOSITE  MATERIALS  INC 

D  SHORTT 
19105  63  AVENE 
PO  BOX  25 

ARLINGTON  WA  98223 

1  COMPOSITE  MATERIALS  INC 

R  HOLLAND 
11  JEWEL  COURT 
ORINDACA  94563 

1  COMPOSITE  MATERIALS  INC 
C  RILEY 

14530  S  ANSON  AVE 
SANTA  FE  SPRINGS  CA  90670 

2  COMPOSIX 
D  BLAKE 

L  DIXON 

120  0  NEEL  DR 

HEBRUN  OHIO  43025 

4  CYTEC  FIBERITE 
R DUNNE 
DKOHU 
MGELIO 
RMAYHEW 
1300  REVOLUTION  ST 
HAVRE  DE  GRACE  MD  21078 

2  SIMULA 

J  COLTMAN 
RHUYETT 
10016  S51ST  ST 
PHOENIX  AZ  85044 

1  SIOUX  MFG 

BKRIEL 
PO  BOX  400 
FT  TOTTEN  ND  58335 
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PROTECTION  MATERIALS  INC 

1 

BATTELLE 

M  MILLER 

CR  HARGREAVES 

FCRnJ.EY 

505  KING  AVE 

14000  NW  58  CT 

COLUMBUS  OH  43201-2681 

MIAMI  LAKES  FL  33014 

2 

BATTELLE  NATICK  OPERATIONS 

FOSTER  MILLER 

J  CONNORS 

JJGASSNER 

BHALPIN 

MROYLANCE 

209  W  CENTRAL  ST 

WZUKAS 

STE302 

195  BEAR  HILL  RD 

NATICK  MA  01760 

WALTHAM  MA  02354-1196 

1 

BATTELLE  NW  DOE  PNNL 

ROM  DEVELOPMENT  CORP 

THALLMSK231 

ROMEARA 

BATTELLE  BLVD 

136  SWINEBURNE  ROW 

RICHLAND  WA  99352 

BRICK  MARKET  PLACE 

NEWPORT  RI 02840 

3 

PACMC  NORTHWEST  LAB 

TEXTRON  SYSTEMS 

M  SMITH 

GVAN  ARSDALE 

T FOLTZ 

R  SHIPPELL 

M  TREASURE 

PO  BOX  999 

201  LOWELL  ST 

RICHLAND  WA  99352 

WILMINGTON  MA  08870-2941 

1 

ARMTEC  DEFENSE  PRODUCTS 

IPS  GLASS 

SDYER 

L CARTER 

85  901  AVE  53 

PO  BOX  260 

PO  BOX  848 

SLATER  RD 

COACHELLA  CA  92236 

SLATER  SC  29683 

2 

ADVANCED  COMPOSITE 

0  GARA  HESS  &  EISENHARDT 

MATES  CORP 

M  GILLESPIE 

PHOOD 

9113LESAINTDR 

J  RHODES 

FAIRFIELD  OH  45014 

1525  S  BUNCOMBE  RD 

MnUKEN  RESEARCH  CORP 

HKUHN 

2 

GREER  SC  29651-9208 

GLCCINC 

M  MACLEOD 

JRAY 

PO  BOX  1926 

M BRADLEY 

SPARTANBURG  SC  29303 

103  TRADE  ZONE  DR 

CONNEAUGHT  INDUSTRIES  INC 

STE  26C 

WEST  COLUMBIA  SC  29170 

J  SANTOS 

PO  BOX  1425 

2 

AMOCO  PERFORMANCE 

COVENTRY  RI  028 16 

PRODUCTS 

MMICHNOJR 

J  BANISAUKAS 

4500  MCGINNIS  FERRY  RD 
ALPHARETTA  GA  30202-3944 
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I  SAIC 

M  PALMER 

2109  AIR  PARK  ROSE 

ALBUQUERQUE  NM  87 106 

1  SAIC 

G  CHRYSSOMALUS 
3800  W  80TH  ST  STE  1090 
BLOOMINGTON  MN  55431 

1  AAI  CORPORATION 

TGSTASTNY 
PO  BOX  126 

HUNT  VALLEY  MD  2 1 030-0 1 26 

1  JOHN  HEBERT 

PO  BOX  1072 

HUNT  VALLEY  MD  21030-0126 

12  ALLIANT  TECHS  YSTEMS  INC 

C  CANDLAND 
C AAKHUS 
R  BECKER 
BSEE 

N VLAHAKUS 

RDOHRN 

SHAGLUND 

DHSHER 

W  WORRELL 

R  COPENHAFER 

MHISSONG 

DKAMDAR 

600  2ND  ST  NE 

HOPKINS  MN  55343-8367 

3  ALLIANT  TECHS  YSTEMS  INC 

J  CONDON 
ELYNAM 
J GERHARD 
WVOl  16  STATE  RT  956 
PO  BOX  210 

ROCKET  CENTER  WV  26726-0210 

1  APPLIED  COMPOSITES 

WGRISCH 
333  NORTH  SIXTH  ST 
ST  CHARLES  IL  60174 


1  PROJECTILE  TECHNOLOGY  INC 

515  GILES  ST 

HAVRE  DE  GRACE  MD  21078 

1  CUSTOM  ANALYTICAL 
ENG  SYS  INC 

A  ALEXANDER 
13000  TENSOR  LNNE 
FUNTSTONE  MD  2 1 530 

2  LORAL  VOUGHT  SYSTEMS 
G  JACKSON 

KCOOK 

1701  W  MARSHALL  DR 
GRAND  PRAIRIE  TX  75051 

5  AEROJET  GEN  CORP 
D  PBLLASCH 
T COULTER 
C  FLYNN 
D  RUBAREZUL 
M  GREINER 

1 100  WEST  HOLLYVALE  ST 
AZUSA  CA  91702-0296 

3  HEXCELINC 
RBOE 
FPOUCELU 
J  POESCH 
PO  BOX  98 
MAGNA  UT  84044 

3  HERCULES  INC 
G  KUEBEIER 
J  VERMEYCHUK 
BMANDERVBLLEJR 
HERCULES  PLAZA 
WILMINGTON  DE  19894 

1  BRIGS  COMPANY 
J  BACKOFEN 
2668  PETERBOROUGH  ST 
HERNDON  VA  22071-2443 

1  ZERNOW  TECHNICAL  SERVICES 
LZERNOW 

425  W  BONITA  AVE  STE  208 
SAN  DIMAS  CA91773 
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2  OUN  CORPORATION 
FUNCHBAUGHDIV 
E  STEINER 
B  STEWART 
PO  BOX  127 
RED  UON  PA  17356 

1  OUN  CORPORATION 

L  WHITMORE 
10101  9TH  ST  NORTH 
ST  PETERSBURG  FL  33702 

1  DOWUT 

STTDRICK 
15STERUNGDR 
WALLINGFORD  CT  06492 

5  SIKORSKY  AIRCRAFT 

GJACARUSO 
TCARSTENSAN 
BKAY 

S  GARBO  MSS330A 
JADELMANN 
6900  MAIN  ST 
PO  BOX  9729 

STRATFORD  CT  06497-9729 

1  PRATT  &  WHITNEY 

DHAMBRICK 
400  MAIN  STMS  114  37 
EAST  HARTFORD  CT  06108 

1  AEROSPACE  CORP 
G  HAWKINS  M4  945 

2350  E  EL  SEGUNDO  BLVD 
EL  SEGUNDOCA  90245 

2  CYTECFIBERnE 
MUN 

WWEB 

1440  N  JCRAEMER  BLVD 
ANAHEIM  CA  92806 

1  HEXCEL 

TBrrZER 

11711  DUBLIN  BLVD 
DUBLIN  CA  94568 


1  BOEING 
RBOHLMANN 

PO  BOX  516  MC  5021322 
ST  LOUIS  MO  63 1 66-05 1 6 

2  BOEING  DEFENSE 
&  SPACE  GRP 

W  HAMMOND 

J  RUSSELL 

S4X55 

PO  BOX  3707 

SEATTLE  WA  98124-2207 

2  BOEING  ROTORCRAFT 

P  MINGURT 
P HANDEL 
800  B  PUTNAM  BLVD 
WALLINGFORD  PA  19086 

1  BOEING 

DOUGLAS  PRODUCTS  DIV 
LJ  HART  SMITH 
3855  LAKEWOOD  BLVD 
D8000019 

LONG  BEACH  CA  90846-0001 

1  LOCKHEED  MARTIN 

S REEVE 
8650  COBB  DR 
D73  62MZ0648 
MARIETTA  GA  30063-0648 

1  LOCKHEED  MARTIN 

SKUNK  WORKS 
D  FORTNEY 
1011  LOCKHEED  WAY 
PALMDALE  CA  93599-2502 

1  LOCKHEED  MARTIN 

RFIEDDS 
1195IRWDJCT 
WINTER  SPRINGS  FL  32708 

1  MATERIALS  SCIENCES  CORP 

B  W  ROSEN 

500  OFHCE  CENTER  DR  STE  250 
FORT  WASHINGTON  PA  19034 
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1  NORTHROP  GRUMMAN  CORP 
ELECTRONIC  SENSORS  & 
SYSTEMS  DIV 
ESCHOCH 
MAILSTOPV  16 

1 745 A  WEST  NURSERY  RD 
UNTHICUMMD  21090 

2  NORTHROP  GRUMMAN 
ENVIRONMENTAL  PROGRAMS 
ROSTCRMAN 

A  YEN 

8900  E  WASHINGTON  BLVD 
PICO  RIVERA  CA  90660 

1  UNITED  DEFENSE  LP 

D  MARTIN 
PO  BOX  359 

SANTA  CLARA  CA  95052 

1  UNITED  DEFENSE  LP 
G  THOMAS 
POBOX  58123 
SANTA  CLARA  CA  95052 

2  UNITED  DEFENSE  LP 
R BARRETT 
VHORVATICH 
MAIL  DROP  M53 
328WBROKAWRD 
SANTA  CLARA  CA  95052-0359 

3  UNITED  DEFENSE  LP 
GROUND  SYSTEMS  DIVISION 
M  PEDRAZZI  MAIL  DROP  N09 
A  LEE  MAIL  DROP  Nil 

M  MACLEAN  MAIL  DROP  N06 
1205  COLEMAN  AVE 
SANTA  CLARA  CA  95052 

4  UNITED  DEFENSE  LP 
4800  EAST  RIVER  RD 
R  BRYNSVOLD 
PJANKEMS170 

T  GIOVANETTI MS236 
B  VANWYKMS389 
MINNEAPOUS  MN  55421-1498 


2  GENERAL  DYNAMICS 
LAND  SYSTEMS 
DREES 
MPASIK 
PO  BOX  2074 
WARREN  MI  48090-2074 

1  GENERAL  DYNAMICS 

LAND  SYSTEMS 
DBARTLE 
POBOX  1901 
WARREN  MI  48090 

1  GENERAL  DYNAMICS 

LAND  SYSTEMS 
MUSKEGON  OPERATIONS 
W  SOMMERS  JR 
76  GETTY  ST 
MUSKEGON  MI  49442 

1  GENERAL  DYNAMICS 
AMPHIBIOUS  SYS 
SURVIVABILITY  LEAD 
G  WALKER 

991  ANNAPOLIS  WAY 
WOODBRIDGE  VA  22191 

5  INST  FOR  ADVANCED  TECH 

TKIEHNE 
HFAIR 
PSULUVAN 
WREINECKE 
IMCNAB 

4030  2  W  BRAKER  LN 
AUSTIN  TX  78759 

2  CIVIL  ENGR  RSCH  FOUNDATION 
H  BERNSTEIN  PRESIDENT 

R BELLE 

1015  15THSTNWSTE600 
WASHINGTON  DC  20005 

1  ARROW  TECH  ASSO 

1233  SHELBURNE  RD  SEE  D  8 
SOUTH  BURLINGTON  VT 
05403-7700 
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1  CONSULTANT 

REICHELBERGER 
409  W  CATHERINE  ST 
BEL  AIR  MD  21014-3613 

1  UCLA  MANE  DEPT  ENGR IV 
H  THOMAS  HAHN 

LOS  ANGELES  CA  90024-1597 

2  UNIV  OF  DAYTON  RESEARCH  INST 
RANYKIM 

AJITKROY 

300  COLLEGE  PARK  AVE 
DAYTON  OH  45469-0168 

1  MTT 

PLAGACE 
77  MASS  AVE 
CAMBRIDGE  MA  01887 

1  nr  RESEARCH  CTR 

DROSE 
201  MILL  ST 
ROME  NY  13440-6916 

1  GEORGIA  TECH  RESEARCH  INST 

GEORGIA  INST  OF  TECHNOLOGY 
PFRIEDERICH 
ATLANTA  GA  30392 

1  MICHIGAN  ST  UNIV 

RAVERILL 
3515  EB  MSMDEPT 
EAST  LANSING  MI  48824- 1 226 

1  UNIV  OF  KENTUCKY 

LPENN 

763  ANDERSON  HALL 
LEXINGTON  KY  40506-0046 

1  UNIV  OF  WYOMING 

D  ADAMS 
PO  BOX  3295 
LARAMIE  WY  82071 


NO.  OF 

COPIES  ORGANIZATION 


1  UNIV  OF  UTAH 

DEPT  OF  MECH  &  INDUSTRIAL  . 
ENGR 

S  SWANSON 

SALT  LAKE  CITY  UT  841 12 

2  PENNSYLVANIA  STATE  UNIV 
RMCNITT 

CBAKIS 

227  HAMMOND  BLDG 
UNIVERSITY  PARK  PA  16802 

1  PENNSYLVANIA  STATE  UNIV 

RENATA  S  ENGEL 
245  HAMMOND  BLDG 
UNIVERSITY  PARK  PA  16801 

1  PURDUE  UNIV 

SCHOOL  OF  AERO  &  ASTRO 
CTSUN 

W  LAFAYETTE  IN  47907-1282 

1  STANFORD  UNIV 

DEPARTMENT  OF  AERONAUTICS 
AND  AEROBALUSTICS 
DURANT  BUILDING 
STSAI 

STANFORD  CA  94305 

1  UNIV  OF  DAYTON 

JMWHTINEY 
COLLEGE  PARK  AVE 
DAYTON  OH  45469-0240 

7  UNIV  OF  DELAWARE 

CTR  FOR  COMPOSIIE  MATRLS 

J  GILLESPIE 

MSANTARE 

GPALMESE 

S  YARLAGADDA 

S ADVANI 

DHEIDER 

DKUKICH 

201  SPENCER  LABORATORY 
NEWARK  DE  19716 
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UNIV  OF  ILLINOIS 
AT  URB  ANA  CHAMPAIGN 
NATL  CTR  FOR  COMPOSITE 
MATERIALS  RESEARCH 
2 1 6  TALBOT  LABORATORY 
J  ECONOMY 
104  S  WRIGHT  ST 
URBANABL  61801 

THE  UNIV  OF  TEXAS 
AT  AUSTIN 

CTR  FOR  ELECTROMECHANICS 

J  PRICE 

A  WALLS 

JKITZMILIER 

10100  BURNET  RD 

AUSTIN  TX  78758-4497 

VA  POLYTECHNICAL  INST 
STATE  UNIV 
.DEPTOFESM 
MWHYER 
KREEFSNIDER 
R  JONES 

BLACKSBURG  VA  24061-0219 

NORTH  CAROLINA  STATE 
UNIV 

CIVIL  ENGINEERING  DEPT 
W  RASDORF 
PO  BOX  7908 
RALEIGH  NC  27696-7908 

UNIV  OF  MARYLAND 

DEPT  OF  AEROSPACE  ENGINEERING 

ANTHONY]  VEZINI 

COLLEGE  PARK  MD  20742 

DREXELUNIV 
ALBERT  SD  WANG 
32ND  AND  CHESTNUT  STREETS 
PHILADELPHIA  PA  19104 

SOUTHWEST  RSCH  INST 
ENGR  &  MATL  SCIENCES  DIV 
JRIEGEL 

6220CULEBRARD 
PO  DRAWER  28510 
SAN  ANTONIO  TX  78228-0510 


ABERDEEN  PROVING  GROUND 

1  COMMANDER 

US  ARMY  MATERIEL  SYS 
ANALYSIS 
P  DIETZ 

392  HOPKINS  RD 
AMXSYTD 
APGMD  21005-5071 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 

AMSRLOPAPL 

APG  MD  21005  5066 

115  DIRUSARL 
AMSRLCI 
AMSRLCIH 
WSTUREK 
AMSRLCIS 
AMARK 
AMSRLCSIOFI 
M  ADAMSON 
AMSRLSLB 
J  SMITH 
AMSRLSLBA 
AMSRL  SL  BL 
DBELY 
R HENRY 
AMSRL  SLBG 
A  YOUNG 
AMSRL SL I 
AMSRL  WMB 
A HORST 
E  SCHMIDT 
AMSRL  WM  BA 
W  DAMICO 
F BRANDON 
AMSRL  WMBC 
PPLOSTINS 
DLYON 
J  NEWEL 
S  WDLKERSON 
A  ZIELINSKI 
AMSRL  WMBD 
B  FORCH 
RFIFER 

R  PESCE  RODRIGUEZ 
BRICE 
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ABERDEEN  PROVING  GROUND  (CONT> 

AMSRLWMBE 

GWKEN 

CLEVERITT 

DKOOKER 

AMSRLWMBR 
C  SHOEMAKER 
JBORNSTEIN 

AMSRLWMM 
D  VDECHNICKI 
GHAGNAUER 
J  MCCAULEY 
B TANNER 

AMSRLWMMA 
RSHUFORD 
PTOUCHET 
N  BECK  TAN 
D  FLANAGAN 
LGHIORSE 
D  HARRIS 
S  MCKNIGHT 
PMOY 
SNGYUEN 
P  PATTERSON 
G  RODRIGUEZ 
ATEETS 
RYIN 

AMSRLWMMB 
BFINK 
J  BENDER 
TBLANAS 
TBOGETn 
RBOSSOU 
L BURTON 
KBOYD 
S  CORNEUSON 
PDEHMER 
R  DOOLEY 
WDRYSDALE 
GGAZONAS 
S  GHIORSE 
D  GRANVILLE 
D  HOPKINS 
CHOPPEL 
D  HENRY 
RKASTE 
MKLUSEWTTZ 
MLEADORE 
RUEB 


ABERDEEN  PROVING  GROUND  (CON'D 

AMSRLWMMB 
ERIGAS 
J  SANDS 
D  SPAGNUOLO 
W  SPURGEON 
JTZENG 
E  WETZEL 
AABRAHAMIAN 
M  BERMAN 
AFRYDMAN 
TU 

W  MCINTOSH 
ESZYMANSKI 
AMRSLWMMC 
J  BEATTY 
JSWAB 
ECHIN 

J  MONTGOMERY 
AWERESCZCAK 
JLASALVIA 
J  WELLS 
AMSRLWMMB 
WROY 
S  WALSH 
AMSRLWMT 
B  BURNS 
AMSRLWMTA 
WGHUCH 
T HAVEL 
J  RUNYEON 
M  BURKINS 
EHORWATH 
B  GOOCH 
W  BRUCHEY 
AMSRLWMTC 
R  COATES 
AMSRLWMTD 
A  DAS  GUPTA 
THADUCH 
TMOYNIHAN 
F  GREGORY 
ARAJENDRAN 
M  RAFTENBERG 
MBOTELER 
TWEERASOORIYA 
D  DANDEKAR 
A  DIETRICH 


43 


NO.  OF 

COPffiS  ORGANIZATION 


ABERDEEN  PROVING  GROUND  TCONT) 

AMSRLWMTE 
ANIILER 
J  POWELL 
AMSRLSS  SD 
H  WALLACE 
AMSRLSS  SER 
R CHASE 
AMSRLSS  SEDS 
RREYZER 
R  ATKINSON 
AMSRLSEL 
RWEINRAUB 
J  DESMOND 
D  WOODBURY 
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1  RMARUN 
MERL 
LTD 

TAMWORTHRD 
HERTFORD  SG13  7DG 
UNITED  KINGDOM 

1  PWLAY 

SMC  SCOTLAND 

DERA  ROSYTH 

ROS YTH  ROYAL  DOCKYARD 

DUNFERMLINE  FIFE  KY  1 1  2XR 

UNITED  KINGDOM 

1  TGOTTESMAN 

CIVIL  AVIATION  ADMINISTRATION 
POBOX8 

BEN  GURION  INTERNE  AIRPORT 
LOD  70150  ISRAEL 

1  S ANDRE 

AEROSPATIALE 
ABTECCRTEMD132 
316  ROUTE  DE  BAYONNE 
TOULOUSE  31060 
FRANCE 

1  J  BAUER 

DAIMLER  BENZ  AEROSPACE 
D  81663  MUNCHEN 
MUNICH 
GERMANY 

3  DRA  FORT  HALSTEAD 

PETER  N  JONES 
DAVID  SCOTT 
MIKE  HINTON 

SEVEN  OAKS  KENT  TN  147BP 
UNITED  KINGDOM 

1  FRANCOIS  LESAGE 

DEFENSE  RESEARCH  ESTAB 
VALC  ARTIER 
PO  BOX  8800 

COURCELETTE  QUEBEC  COA 
IRO  CANADA 


2  ROYAL  MCLITARY  COLLEGE  OF 

SCIENCE  SHRIVENHAM 
DBULMAN 
B  LAWTON 

SWINDON  WILTS  SN6  8LA 
UNITED  KINGDOM 

1  SWISS  FEDERAL  ARMAMENTS 

WKS 

WALTER  LANZ 
ALLMENDSTRASSE  86 
3602  THUN 
SWITZERLAND 

1  PROFESSOR  SOL  BODNER 

ISRAEL  INST  OF 
TECHNOLOGY 

FACULTY  OF  MECHANICAL  ENGR 
HAIFA  3200  ISRAEL 

1  DSTO  MATERIALS  RSRCH  LAB 

DR  NORBERT  BURMAN  NAVAL 
PLATFORM  VULNERABILITY  SHIP 
STRUCTURES  &  MATERIALS  DIV 
PO  BOX  50 

ASCOT  VALE  VICTORIA 
AUSTRALIA  3032 

1  PROFESSOR  EDWARD  CELENS 

ECOLE  ROYAL  MUJTAIRE 
AVE  DE  LA  RENAISSANCE  30 
1040  BRUXELLE 
BELGIQUE 

1  DEF  RES  ESTABUSHMENT 

VALCARTIER 
ALAIN  DUPUIS 

2459  BOULEVARD  PIE  XI  NORTH 
VALCARTIER  QUEBEC 
CANADA 

PO  BOX  8800  COURCELETTE 
GOA  IRO  QUEBEC  CANADA 
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INSnrUT  FRANCO  ALLEMAND  DE 

RECHERCHES  DE  SANIT  LOUIS 

DEMARCGIRAUD 

RUE  DU  GENERAL  CASSAGNOU 

BOITE  POSTALE  34 

F  68301  SAINT  LOUIS  CEDEX 

FRANCE 

JMANSON 
ECOLE  POLYTECH 
DMXLTC 

CH  1015  LAUSANNE  SWITZERLAND 

TNO  PRINS  MAURTTS  LAB 
ROB  USSELSTEIN 
LANGE  KLEIWEG  137 
POBOX45 
2280  AA  RUSWUK 
THE  NETHERLANDS 

FOA  NATL  DEFENSE 
RESEARCH  ESTAB 
BOJANZON 
RHOLMUN 

DIR  DEPT  OF  WEAPONS  & 
PROTECTION 
S  172  90  STOCKHOLM 
SWEDEN 

DEFENSE  TECH  &  PROC  AGENCY 
GRND 

ICREWTHER 
GENERAL  HERZOG  HAUS 
3602  THUN 
SWITZERLAND 

MINISTRY  OF  DEFENCE 
RAFAEL 

MEIRMAYSELESS 
ARMAMENT  DEVELOPMENT  AUTH 
PO  BOX  2250 
HAIFA  31021  ISRAEL 

AKE  PERSSON 
DYNAMEC  RESEARCH  AB 
PARADISGRND  7 
S  151  36  SODERTAUE 
SWEDEN 


1  ERNST  MACH  ESISTITUT  EMI 

DIRECTOR 
HAUPTSTRASSE  18 
79576  WELL  AM  RHEIN 
GERMANY 

1  ERNST  MACH  INSHTUT  EMI 

ALOIS  STTLP 
ECKERSTRASSE  4 
7800  FREIBURG 
GERMANY 

1  IRHANSPASMAN 

TNO  DEFENSE  RESEARCH 
POSTBUS  6006 
2600  JA  DELFT 
THE  NETHERLANDS 

1  BITAN  HIRSCH 

TACHKEMONY  ST  6 

NETAMUA42611 

ISRAEL 

I  MANFRED  HELD 

DEUTSCHE  AEROSPACE  AG 
DYNAMICS  SYSTEMS 
PO  BOX  1340 

D  86523  SCHROBENHAUSEN 
GERMANY 
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